Quantum description and properties of electrons emitted from pulsed nanotip electron 

sources 
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We present a quantum calculation of the electron degeneracy for electron sources. We explore 
quantum interference of electrons in the temporal and spatial domain and demonstrate how it can 
be utilized to characterize a pulsed electron source. We estimate effects of Coulomb repulsion on 
two-electron interference and show that currently available nano tip pulsed electron sources operate 
in the regime where the quantum nature of electrons can be made dominant. 
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I. INTRODUCTION 

Ever shorter synchronous laser and electron pulses are 
being developed [ll-ll| as probes to investigate the spatial 
and temporal behavior of atoms, molecules and nanos- 
tructures. Although attosecond physics is already an 
exciting and active area of research, no electron pulses 
shorter than ten femtoseconds that are deliverable on a 
target of choice have been realized yet. After the initial 
generation of electrons, the short pulses rapidly disperse 
due to their energy spread. The energy spread is inherent 
to the Heisenberg uncertainty principle and pulse com- 
pression is thus needed to deliver short pulses on a target. 
Pulse compression techniques such as RF pill boxes 
and shaped laser modes [l| are being explored for this 
purpose. 

In this paper ultrafast electron dynamics is considered. 
For one-electron pulses, the pulse duration is limited 
by the uncertainty principle. This is expected to pose 
no principle obstacle for current technologies to reach 
deep into the attosecond regime [l|. For multi-electron 
pulses, wave packet compression is ultimately limited by 
Coulomb repulsion and the Pauli Exclusion Principle. To 
address the question how dynamics of wave packets is in- 
fluenced by the Exclusion Principle, a theoretical frame- 
work is developed in this paper. Using this approach, es- 
timates of the effects of electron degeneracy are made and 
compared to estimates of the effects of Coulomb repul- 
sion. We show that at the source degeneracies of about 
10~^ are expected in current pulsed electron sources ig- 
noring Coulomb repulsion. The effects of Coulomb re- 
pulsion are comparable in magnitude to the dynamical 
effects of degeneracy. The degeneracy of 10^^ is three 
orders of magnitude higher than previously reported for 
any free electron experiment with a continuous electron 
source. This is in spite of the fact that continuous field 
emission tips are considered to be very bright sources and 
a workhorse for electron microscopy Pulsed sources 
are thus expected to increase the instantaneous bright- 
ness and degeneracy considerably. Electron pulse com- 
pression may increase the degeneracy even further, and 
it is expected that the Pauli Exclusion Principle needs to 



be accounted for in the dynamics of ultrafast free electron 
pulse studies. 

Our theoretical approach is designed to handle the evo- 
lution of electron wave packets that evolve continuously 
in time from the initial to the final state. Rather then 
studying the process of electron emission from a nano 
tip, we assume that a certain multi electron state was 
emitted by the source. We then are interested in its 
time evolution, taking into account electron dispersion 
and Coulomb repulsion. The approach leads to expres- 
sions that are related to the intensity correlation func- 
tions that one finds for the Hanbury-Brown Twiss ef- 
fect j7|] , but provides a time-dependent picture as well as 
the usual spatial picture. An example is discussed of the 
first order electron correlation signal measured with one 
detector that can be used to determine the magnitude of 
the electron degeneracy. That this is possible relies on 
the consideration of a pulsed source as opposed to the 
usual continuous source. Continuous sources were used 
in the historic optical experiment of Hanbury-Brown and 
Twiss , and in the more recent experiment that demon- 
strated the same effect for free electrons @. 

This paper is organized as follows. In Section |TT] the 
description formalism is introduced. A model for the den- 
sity matrix of the source is discussed, and the degener- 
acy is defined. The values of the degeneracy for a pulsed 
electron source are calculated in Section IIII Bl using ex- 
perimentally feasible parameters. In Section lTVl auantum 
interference effects for two-electron states are discussed. 
The effects of Coulomb repulsion are estimated in Sec- 
tion |Vl Summary, conclusions and outlook are provided 
in Section IVTl 
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II. DESCRIPTION FORMALISM 

Electron sources to be considered throughout this pa- 
per can be sketched as consisting of a metallic nano tip 
irradiated by a pulsed laser. The duration of each laser 
pulse is relatively short - on the order of 50 femtoseconds. 
As a result, the tip is exposed to very strong electric fields 
over a very short period of time. The intensity of the laser 
can be adjusted and thus the average number of photo- 
electrons emitted by the tip can be controlled ranging 
from just a single electron to hundreds of electrons per 
pulse 0. 
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We are interested in the scenario when, on average, 
only a few electrons are emitted per pulse. This would 
be the regime where the quantum nature of electrons is 
dominant and where a quantum theory describing the 
source's statistics is needed. 



A. Second Quantization 

The standard formalism of quantum mechanics applies 
to systems with fixed number of particles. However, this 
is not the case when one is interested in electron emis- 
sion sources. The number of electrons emitted by the 
tungsten nano tip fluctuates from pulse to pulse. In this 
situation, the second quantization formalism is a natural 
way to account for the fluctuations within the quantum 
paradigm. 

Let us briefly review the basics of the second quantized 
description. Instead of keeping track of individual elec- 
trons we will use the notion of electron field. The field 
is represented by the field operators ^'^(x), *I'(x) defined 
continuously in all points of space x. They obey the 
following anticommutation relations {'I'(x), ^''^(x')} = 
(5(x — x') - a direct consequence of the fermionic nature 
of electrons. It is also convenient to define creation and 
annihilation operators cj^ and Ck that are connected to 
the field operators via the Fourier transform: 

^'(x) = ^ / dke^'^-^Ck 

yV J-oo 

= / dke-'^-ct, (1) 

here is a quantization volume. Like the field operators, 
the creation and annihilation operators also anticommute 
i.e. {ck,c[,} = <5(k-k'). 

The next step is to introduce the Hilbert space that 
encompasses all possible multi-electron states. It is also 
known as the Fock space and it is constructed as a direct 
sum of tensor products of single-particle Hilbert spaces. 
Among all of the Fock states, the most basic one is the 
vacuum |0). It corresponds to a situation when no elec- 
trons are emitted by the tip. Acting on the vacuum with 
the creation operator cj^ will generate a particle (an elec- 
tron) in a state described by the index k ^1^. In other 
words, cj^.|0) — |0 . . . OlkO . . • 0) = |lk), that is, an elec- 
tron in the state |k) is created from the vacuum. On 
the other hand, by acting with the annihilation opera- 
tor Ck on the one-electron Fock state | Ik) a particle will 
be destroyed and the resulting state will be the vacuum 
again. Analogously, we interpret the action of the field 
operator 'I''f(x) upon the vacuum state |0) as a creation 
of a particle in the position x, or equivalently an electron 
in the state |x), i.e. ^'''(x)|0) = |lx). Lastly, because 
the Pauli exclusion principle for electrons imposes that 
{C]{^,C]{^,} — no two electrons can be found in a one- 
particle state described by the same state index and thus 



states like |2k) are forbidden. However, the states like 
|lk • • • Ik') are allowed. 

B. State Description: Density Matrices 

Having previously defined electron field operators and 
the Fock space they act upon we now need to assign a 
quantum state to our electron source in order to be able 
to make quantitative predictions about electron statis- 
tics. Let us denote the total number of laser pulses that 
strike the tip of the source by M . Suppose, that after 
the z-th laser pulse a bunch of electrons in the state 
is emitted. Assuming that Nmax is the maximal number 
of electrons that can be emitted per pulse, one can write 
down l^*;) as follows, 

\^,) ^ aoM + ai.^ll'^') + . . . + «jV™...Ka,), (2) 

where |0) is the vacuum state (no electrons emmited), 
If^') is a state with one electron, [i V"^' ) is a state that 
contains precisely Nmax electrons [i2[, and \aj^i\'^,j = 
0, . . . , Nmax is a probability to detect j electrons after 
the i-th pulse. 

We are now going to make a major assumption about 
the nature of our source. We assume that the proba- 
bilities to detect 0, 1, ... , Nmax electrons do not change 
from pulse to pulse i.e. is the same for all i for a 

given value of j. The assumption can be backed up ex- 
perimentally. Indeed, at currently achievable laser rep- 
etition rates of several MHz, the time interval between 
two successive electron emission events is several orders 
of magnitude larger than the time scale of relevant dy- 
namics within the emitted multi-electron states [1, |Tlj . 
In other words, the states should not be affected by 
the history of prior electron emission events. Therefore, 
provided that all experimental parameters relevant to the 
emission process are stabilized, our assumption should be 
reasonable. 

In general, however, the coefficients aj,i are complex 
numbers. This means that in order to fully control the 
states I'i'i) one must not only control the absolute values 
of aj^i but also their phases. The later is extremely chal- 
lenging; such degree of control is not present in current 
experiments. Thus, we have to assume that aj^i acquire 
essentially random phases from pulse to pulse. 

Because the states will all have different relative 
phases for each laser pulse i, no pure state can be assigned 
to our source. Instead the state of the source is described 
by a density matrix p^, 

1 

= mEi*«)(*'^I- (3) 

1=1 

Expanding this sum using the definition of j^I'i) in Eq.([5]) 
we obtain, 

Ps — pjPj + cross terms, (4) 

3=0 
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here pj = \ctj\'^,j = 0,Njnax is the probabihty to emit 
exactly j electrons (note that by our assumption these 
probabilities are the same for every pulse i so we have 
omitted the index i all together), pj — is a j- 

electron density matrix and its detailed structure will be 
discussed shortly. The cross terms in Eq.Q are of the 
form p„in — Ito*^') (n"^'!, m, n = 0, ■ ■ ■ , N^ax where \m'^'-) 
and are TO-(n-) electron states respectively. Each 

cross term pmn comes with a pre factor 



y/PmP 



M 



1=1 



(5) 



where PmiPn) are m-(n-)electron generation probabil- 
ities, 4'i{m,rL) is the random relative phase produced 
after the ?-th laser pulse. Assuming that the phases 
4>i{m,n) are distributed uniformly over [0, 27r] interval 
and M — cx) we conclude, using the central limit the- 
orem, that Pmn — ^ 0. In other words, the quantum co- 
herence terms between various electron states vanish due 
to our inability to control their relative phases between 
pulses. 

Taking the above discussion into account we can write 
the density matrix ps as 



N„ 



E 



PjPj. 



(6) 



Let us now discuss the structure of j-electron states in 
more details. The simplest case is j = 0, i.e. the vacuum. 
In this case po = |0)(0|. Let us further assume, since 
Nmax can be controlled experimentally e.g. by tuning the 
intensity of the laser, that Nmax = 2. This will simplify 
our calculations on one hand, yet, on the other, it will 
still provide us with an insight into the multi-electron 
quantum phenomenon. 



C. One-electron States 

Formally, we can write the one-electron density matrix 
pi = But how does one model Imag- 

ine for a second that we have chosen = |lk). In 
other words, wc have assigned an eigenstate |k) of the 
one-electron momentum operator to our particle. If we 
now measure the momentum of the electron we should 
get the value of k = (/cj,, ky,kz) as an outcome. Also, the 
uncertainty {Ak^, Aky, Akz) of the measurement should 
be zero. On the other hand, the emission tip occupies a 
finite volume AV = AxAyAz. This provides limitations 
on the degree of spatial localization of the electron. Thus, 
using the Heisenberg uncertainty principle, we conclude 
that the uncertainty in measuring the momentum of the 
electron Ak^AkyAk^ > We immediately see that 

our initial state assignment jl*^') = |lk) contradicts the 
uncertainty condition. Therefore, we should rather treat 



the emitted electron as a wave packet i.e. as a superpo- 
sition of different momentum eigenstates: 



+c>o 



dkC(k)|lk), 



(7) 



where |C(k)p is a probability density function that de- 
scribes particular details of the momentum distribution. 
Experimental data ^Ts'] suggest that |C(k)p resembles a 
normal distribution with the mean ko. 



C(k) 



1 



(8) 



where E ^ = 4Afc2 diag[l, 1, 1] and we assume that the 
variances Ak^, Aky, Akz are the same and equal to Ak. 



D. Two-electron States 

In quantum theory identical particles are assumed to 
be indistinguishable. This imposes a symmetry, under 
the operation of permutation, onto systems of two or 
more electrons. More specifically, the indistinguishability 
of electrons requires a multi-electron wave function to be 
antisymmetric under particle permutations. This leads 
to the Pauli exclusion principle (PEP) that forbids two 
electrons to be in the same one particle state. Therefore, 
the PEP should be take into account when constructing 
two-electron states. 

So far we have been ignoring spin degrees of freedom 
that electrons inherently have, concentrating only on the 
spatial part of electron states. However, the total state 
of an electron |4't) is a combination of the spatial and 
spin states, i.e. |*t) = IV'spatmi) I'^spm). The PEP dic- 
tates that for two(or more) electrons the total electron 
state I^'t) must be antisymmetric under particle permu- 
tations. Due to the direct product structure of \^t) only 
two possibilities are present: either \ipspatiai) is antisym- 
metric and \(l)spin) is symmetric or the other way around. 
We will talk about symmetric and antisymmetric spatial 
states in more details in a moment. First, let us concen- 
trate on the spin part of |\E't)- 

It is well known that, there are three symmetric and 
one antisymmetric(think of the triplet and singlet states) 
spin states for two-electron systems. So, in general, |5't) 
can be decomposed into 

Ivf^) ^ iv^^S) « J2 ) + C4|V^) ® (9) 

i=l 

where |0f) denote the spin(symmetric) triplet states, 
10"^'^) is the spin singlet state, and we do not specify the 
explicit structure of the symmetric(antisymmetric) spa- 
tial components |V''^)(|'/''^'^)). In reality, no pure state 
like in Eq. ^ can be assigned to a nanotip electron emis- 
sion source, because experimentally there is no control 
of spin polarization of emitted electrons. Assuming that 
all 4 spin states have equal generation probability(a fully 
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unpolarized spin state), one can write down a density 
matrix describing the total two-electron spatial plus spin 
state: 

3 

PT = \pAS ® J2 I'^f )(<^f I + ® (10) 

i=l 

where pAS = |V'^'^)(V''^'^I and ps = are the 

spatial components of the total density matrix px- More- 
over, since the spin degrees of freedom are neither mea- 
sured nor controlled experimentally one can trace them 
out when describes experimental observations. This leads 
to the following two-electron density matrix, 

3 1 

P2^Tr,pin{pT) ^ ^PAS + ^PS- (11) 

Now let us discuss the choice of ps and pAS ■ We intro- 
duce two one-electron states and \(f>) in the form of 
wave packets given in Eq.Q. We assume that the wave 
packets are centered around ko and kg (ko 7^ kj,) in mo- 
mentum space and have identical variances Ak. Then de- 
fine spatially symmetric and antisymmetric two-electron 
states |2g) and as follows. 



|2|') 



-AS/ 



VW+MW) 
1 



{\ij)<s>\cf,) + \<j>)^m,{u) 

(IV) ® |<^) - 10) ® |V^)),(13) 



where the overlap between the states and {(f)) is given 
by, 



(14) 



Substituting definitions of the states {ip) and |(?!')(Eqs.([7]- 
H])) into Eg. (11211131) we arrive at. 



2sM5)=^±- // '^ki(ik2C±(ki,k2)|lk,,lk.), (15) 



here we have introduced the functions 

kp-kil^ lt^-lt2|'' kp-kal^ k^-kil^ 



(16) 



(27rAfc2)t 

where C_|_ is a symmetric function(C+(ki, k2) 
C_(_(k2,ki)) corresponding to the state |2g), C_ is an an- 
tisymmetric function(C_(ki, k2) = — C_(k2,ki)) corre- 
sponding to the state |2^'g), and. 



N± 



(17) 



The two-electron density matrix can be finally defined by 



E. Electron Statistics 

Electron fields can be interrogated directly, using a 
particle detector that counts the number of electrons in a 
detection volume. A single electron detector provides in- 
formation about electron densities at the location of the 
detector. One can also use multiple detectors in different 
locations to study temporal and spatial correlations of 
electron densities. These data, in principle, is enough to 
gain a complete knowledge of the electron field of inter- 
est. 

The particle density operators at the position x is de- 
fined as /(x) = 'I''l'(x)4'(x). Provided that we know the 
state of the field ps, we can calculate the number of elec- 
trons at x as fT.(x) = Tr(/(x)ps). In a similar fashion the 
number of electrons in a volume AV^ centered around Xq 
can be calculated as, 



n{AV)=Tr{ps (f x¥ {x)} . 



(18) 



Ay 



We recall that an important characteristics of quantum 
particle sources (electron, photon, etc) called the degen- 
eracy [13] is defined as a number of particles per phase- 
space-cell volume. The volume of a phase-space cell is de- 
fined in terms of the Heisenberg uncertainty AxAk = 1. 
Thus, using Eq.(|18|). the electron source degeneracy in 
the vicinity of a point xq reads. 



Tr{ps / d3a;#t(x)*(x)}. 



(19) 



Ax 



substituting the states |2J^g) and |2|') into Eq. pTj) . 



where Ax = -^k is the position uncertainty calculated at 
the source(xo = 0). As a consequence of this definition, 
the electron degeneracy is the largest at the source and 
drops off as one moves away from it. This is true for 
both pulsed and continuous electron sources and is due 
to the dispersion of electron wave packets as they travel 
in space. Indeed, the largest contribution to the degen- 
eracy at Xq, according to Eq. p^ . comes from electron 
wave packets that are centered at xq. But their position 
uncertainty at Xq will be greater than it was at the source 
(Ax) . The number of electrons in the volume Ax at xq 
will appear to be smaller than it was in the same volume 
at the source and, thus, the degeneracy will be less than 
it was at the source. 

The degeneracy attains its largest value when there 
is exactly one electron per phase-space cell (remember, 
electrons are fermions!). This would correspond to the 
degeneracy of 1. In this regime the quantum nature of 
electrons is the most pronounced. In reality, the best re- 
ported value of the degeneracy to date was achieved using 
a continuous electron emission source and is on the order 
of 10^'' [1]. This value is measured at the source and 
even though it is very small it allowed the experimen- 
tal demonstration of the Hanbury-Brown Twiss effect for 
electrons in vacuum. It has also been argued that with 
the degeneracy values of 10~* one could directly observe 
quantum behavior of electrons such as antibunching ■ 
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III. RESULTS: PULSED ELECTRON SOURCE 
DEGENERACY 



and thus easier to interpret. The one-dimensional one- 
particle component of the degeneracy Si reads, 



Calculating the electron degeneracy is now a straight- 
forward task. We just need to combine Eq.® with 
Eg. dlQp . Since the vacuum component of ps does not 
contribute to the degeneracy, we obtain. 



S PiSi +P2S2 



(20) 



where pi,p2 are one-(two-)electron state probabilities, 
and 61,62 are correspondent degeneracy contributions. 
Note that the two-electron degeneracy 62 contains con- 
tributions from spatially symmetric and antisymmetric 
two-electron states according to Ea. ([TT|) . 

So far we have been implicitly assuming that our elec- 
tron source is stationary. In other words, electron statis- 
tics of the source, including the degeneracy, is invariant 
under time translation. This would apply well to elec- 
tron sources operating in a continuous regime. However, 
the source we are interested in here is pulsed and the 
time interval between two consecutive electron emissions 
is several orders of magnitude larger than the duration 
of each pulse. In this situation, the value of the degen- 
eracy will depend not only on the position in space but 
also time. Thus it is imperative that we account for time 
evolution when we calculate the degeneracy. 



Degeneracy Calculation for Non-interacting 
Particles 



To include time evolution into our consideration, we 
need to either specify the time dependence of the electron 
field operator 4'^(x, i), ^(x, t) and keep the state ps in 
Eq. © time independent or keep the field operators time 
independent and instead use the time dependent density 
matrix Ps(t). Here we choose the former approach. 

First, consider non-interacting electrons [19]. The 
Hamiltonian describing the source is then given by. 



n = 



+ 00 



cike(k)cj^Ck, 



(21) 



where e(k) — -1;^. Define the time-dependent electron 
field operators as follows. 



^(x,t) = e"^^*\E'(x,0)e^^* 

f +00 



1 p+00 

W j-00 



1 /■+OC' 

^'t(x,t) ^ — dke-^'^-^e^'C^^'^. (22) 

yV J -00 



One last simplification is made before we calculate the 
time-dependent version of the degeneracy. Let us do the 
calculations for a one-dimensional case first. The general- 
ization to three dimensions is then straightforward. How- 
ever, a one-dimensional solution is algebraically tractable 



Tr{\U){U\ / dx¥{x,t)i!{x,t)} 

J Ax 



1 



hkt.. 



-{erfi^Mt)Ak{x2 )) 

2 TOe 

Hkt^ 



-erfi^M^Ak{xi ))}, (23) 



where erf is the error function, xi = xq — 



Ax 
2 



and 



X2 = xo + ^(here xq and Ax denote the location of 

the detector and its volume), and a(t) = „i2^^2^f^if2 
defines the spatial expansion of a Gaussian wave packet 



in time(Ax(t) 



Aa:(0) 



) jl6|. Note that the largest value 



of the degeneracy is given by srf{^) < 1 and it is at- 
tained when t = 0. This is because the one-electron state 
I lei) at t = is a Gaussian wave packet with the position 
uncertainty Aa; = Therefore, only a part of the wave 
packet will be detected in the volume Aa; around x = 0. 

The result can be easily generalized to three dimen- 
sions, provided that there are no cross correlations be- 
tween kx,ky,kz. In this case the total one-particle con- 
tribution wiU be of the form 6l^{t) = 6fit)6f{t)6f{t), 
where 6\{t),i = x,y,z is given in Ea. (P51) . 

In a similar fashion, using the definition of the two- 
electron density matrix in Ea. ([TT|) . we can calculate the 
symmetric and antisymmetric components of the two- 
electron degeneracy in one dimension. If we define 

U^,t,k) = erf{^Mt)Ak{x-—)), 

me 

6(x,i,fci,fc2) = er/{VMt)Afc(a;-^fci-^) 

Znie 



2\/2Afc 



(fci - fc2)}. 



Ux,tMM) = er/{VMi)Afc(a;-:^^%i^) 



2me 



Va(0 
2\/2Afc 



(fci - fc2)}. 



(24) 



then the two-particle contributions to the degeneracy 
read. 



1 



02iS/AS).^, 2(1 ±(^10)2, 

+ Ci(a;2,i, fco) -6(a;i,i,fco) 

± [6(2^2, i,fco,fco) - £,2{xi,t,ko,k'Q) 

+ 02{x2,t,kn,k'Q) ~ 92{xi,t,ko,k'a)]}, (25) 

where the wave packet overlap {'ip\(j)) is defined in Eg. p^ 
and the plus(minus) sign corresponds to the contribu- 
tion 62s (t) {62 AS (t)) from the symmetric (antisymmetric) 
state \2g){\2'fg)) in Eq. p^ . Let us examine the expres- 
sion for 62{t) carefully. We notice that the difference 
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i,i{x2,t, kg) — ^i{xi,t, ko) = 25i{t, ko) can be interpreted 
as the " one-particle" degeneracy 6i for the electron with 
the mean momentum ?ifco (compare to Eg. (|23|) ) . Simi- 
larly, the difference ^i(a;2, t, fcg) — /cq) — 2Si{t,kQ) 
is the one-electron degeneracy for the electron the mean 
momentum hk'Q. The term in the square brackets in 
Eq. (1^51) contributes to the degeneracy only if the overlap 
between two electron wave packets is non zero. We call 
it the interference term. If the electrons are far apart 
in the momentum space, i.e. \hko — fikgl ^ 1 and Ak 
then they can occupy simultaneously the same positions 
in space and behave like classical particles. In this case 
the two-electron degeneracy reduces to a sum of one par- 
ticle degeneracies and no quantum interference occurs. 
Putting all this together we rewrite Eq. ([25l) as follows, 

S2is/AS){t) = (i±(^|0)2) {'^i(^'feo) + '^ift'feo)(26) 

where we have defined 

2Se{t) = e2{x2,t,ko,k'a)-e2{xi,t,ko,k'o) (27) 
2S^{t) = Ux2,t,ko,k'„) - ^2ixi,t,ko,k'„). (28) 

The generalization of Eq. ipS)) to three dimension is also 
straightforward, 

Slfs/As)it) - j^^^;i;^^mt,koJSl{t,koJSt{t,koJ 
+ 6nt,k'oJSy{t,k',JSt{t,k',J (29) 

± im'remmit) + si{t)si{t)si{t)]}. 

Finally, the total two-electron contribution to the degen- 
eracy for our pulsed electron source reads, 

Sr{t)^\-Sli{t) + l-5lUt). (30) 
B. Results 

Let us now calculate the actual value of the degeneracy 
using realistic numbers for the source parameters. For 
that the following information is needed: 

• the position of the detector with respect to the 
source 

• the direction and the magnitude of the mean elec- 
tron momentum (for both one- and two-electron 
states) 

• the momentum uncertainty in x,y, z directions . 

We assume that the center of the detector is positioned 
on the X axis 0.1m away from the source. We also as- 
sume that the mean electron momenta for one- and two 
electron states have only non-zero components along the 
X axis, i.e. fco^^ = 0, feg, — and /cq^ = 0, fcg^ = 0. The 




8.35 X 10 "'8.4 X 10"' 8.45 x 10"' 8.5 x 10"' 
t,(s) 

FIG. 1: The one-electron degeneracy as a function of time. 
The position of the detector in meters is given by f = 
(0.1,0,0). The mean value of the electron energy is taken 
to be Eo = 400ey {ko, = 1.024 • 10"m~\ fco„ = 0, ko, = 0) 
and the energy spread AE = leV{Ak = 1.28 • lO'^m"^ is 
assumed to be the ScLiiiG for kx 5 ky , ) ■ 

experimental data [l3| suggests that the mean elec- 
tron energy for a pulsed nanotip source Eq can be, for 
example, 400ey and AE sa leV. These translate into 
ko, = 1.024-10"to-i andAk^^y^z = 1.28-lO^m-i. Using 
these values we plot the one-electron, two-electron, and 
total electron 3D degeneracy of the source as a function 
of time in Figs.((l]|4l). Also, to compare the degeneracy 
of our source to the best available continuous sources we 
compute the value of the degeneracy at the source. For 
the above choice of parameters « 0.2, i.e. three or- 
ders of magnitude larger than the value reported in 
Larger values of the degeneracy for pulsed sources are 
expected and can be explained by much higher electron 
current densities generated by pulsed laser light. 

A discussion of the results is in order. First, let us ex- 
amine the one-electron contribution to the degeneracy in 
Fig.©. As one can see the shape of the curve resembles a 
Gaussian. It is due to the fact that our one-electron state 
is a Gaussian wave packet traveling through the detection 
volume. When the center of the one-electron wave packet 
reaches the middle point of the detector(f = (0.1,0,0)) 
the degeneracy attains its maximum value(~ 1.55-10"^'*). 
The reason why this value is ^ 1 is because the wave 
packet spreads out in all spatial directions as a result of 
the propagation, however, our detection volume equals 
the initial spatial spread of the wave packet at t = 0. 
Thus the value of the degeneracy depends on how fast 
the initial wave packet disperses and how far from the 
source the detector is. 

As we have already discussed, the two-electron contri- 
bution to the degeneracy has two components stemming 
from the symmetrized and antisymmetrized states (see 
Fig-©)- The symmetric contribution resembles a gaus- 
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8.4 X 10"' 8.5 X 10 

t,(s) 



FIG. 2: The symmetric (dashed hne) and antisymmetric(soHd 
line) components of the two-electron degeneracy as a func- 
tion of time EQ. (|29p . The position of the detector in me- 
ters is given by r = (0.1,0,0). The x,y,z components of 
the mean momentum vector of the first electron are ko^ — 
1.024 ■ lO^^m"^, fcoj^ = 0, ko^ = and the momentum spread 
Ak = 1.28 ■ 10*m~^ is the same in all directions. For the 
second electron k'g^ = ko^ — 0.5 ■ Ak, k'o^ = 0, fco, = and we 
also assume the same momentum spread Afc = 1.28 ■ 10*m~^ 
in X, y, and z direction. 



sian, meaning that both electrons can be detected at the 
center of the detector at the same time. This is not sur- 
prising since the underlying two-electron wave packet is 
symmetric, i.e. electrons behave like bosons despite the 
fact that they are fermions. On the other hand, the anti- 
symmetric contribution to the degeneracy has two max- 
ima that correspond to the passage of either one of the 
electrons through the detector. In between those max- 
ima there a minimum that is attained when the center 
of the antisymmetrized two-electron wave packet passes 
through the center point of the detector. Why is that? 
Recall that for the two-electron states we have set the 
value of fcg = ko^ — 0.5 • Afc. This means that the an- 
tisymmetric two-electron state consists of two, strongly 
overlapping(in momentum space), wave packets which, 
according to the PEP, may not overlap spatially. In 
other words, the probability to detect two electrons in 
the center of the wave packet drops down to zero. As the 
wave packet propagates the probability is still zero in 
the center but the wings of the wave packet will disperse 
spatially much like in the one-electron case. Thus the 
value of the antisymmetric component of the degeneracy 
will experience a reduction when the center of the two- 
electron packet is at the central point of the detection 
volume. 

The total two-electron degeneracy contribution is a 
weighted sum of the symmetric and antisymmetric com- 
ponents according to Eg. ([50)1 . We plot it on Fig. ([3]) as 




8.4 X 10" 8.5 X 10 

t,(s) 



FIG. 3: The total two-electron degeneracy as a function of 
time defined in EQ. (|30p . The position of the detector, and 
all other relevant parameters are provided in the caption of 
Fig.©. 



l.xlO"'^- 




8.4 X 10"' 8.5 X 10 

t,(s) 



FIG. 4: The total (combined) degeneracy of the electron 
source as a function of time Ea. (|20p . The values of the rele- 
vant parameters are the same as in Fig.([T]) and Fig.©. We 
set the one-electron probability pi — 0.5, and the two-electron 
probability p2 =0.1. 



a function of time. Its distinct feature is a dip in the 
center of the gaussian that is attributed to the destruc- 
tive interference of two electrons constituting the anti- 
symmetric state |2^'g). To observe this two-electron dip 
experimentally one needs a detector with time resolution 
on the order of 20 picoseconds which is a typical time 
resolution of top of the line detectors. 

The total electron degeneracy is a weighted sum of the 
one- and two-electron contributions. It is depicted in 
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Fig. ([4]). The weighting coefficients, i.e. the probabihties 
of one- and two-electron states pi and P2, are chosen 
such that the the one-electron events are dominant i.e. 
Pi = 0.5, p2 — 0.1. Therefore, it is not a surprise that the 
total degeneracy is close in its shape to the one-electron 
component on Fig(IT]). This particular choice of pi and p2 
also implies a high probability of vacuum generation since 
Po = 1 — Pi — P2- In general, pi and p2 can be controlled 
experimentally. For instance, by adjusting the power of 
laser pulses. Then, assuming that electron statistics is 
governed by the Poisson distribution, the two-electron 
events can be made dominant (j)2 > Pi) and one could 
study electron-electron interference by just measuring the 
total degeneracy as a function of time. 

One final remark on how the degeneracy can be mea- 
sured. Of course, in an experiment it is rather the bright- 
ness, defined as the number of particles per unit area per 
unit solid angle per second, that is measured routinely. 
Note, however, that the brightness is closely related to 
the degeneracy [13], and, therefore the value of the de- 
generacy at the source can be determined by measuring 
the brightness. 



IV. QUANTUM INTERFERENCE EFFECTS 
WITH PULSED ELECTRONS 

Quantum interference between two electrons is present 
in the time domain. The consequence of this for degen- 
eracy in the two-electron contribution was discussed in 
Section IIIIBI There, the fact that electrons obey the 
PEP and, therefore, may not be present at the same lo- 
cation simultaneously led to a dip in the degeneracy as a 
function of time(see Fig. ([3])). The same effect leads to a 
spatial interference pattern for two electrons. To observe 
it, one can use two detectors and record the coincidence 
rate of the detectors as a function of their relative dis- 
tance. The resulting spatial pattern will have a dip at the 
point where the relative distance is zero. Note that this 
measurement requires two detectors, whereas to observe 
the dip in the two-electron degeneracy only one detector 
is needed. In this section we discuss two-electron spa- 
tial interference effects that can be observed with one 
and two detectors and study how useful those can be for 
discriminating between states with different number of 
electrons. 



A. Discrimination Between One- and Two-electron 
Components 

The degeneracy of an electron source can be measured 
directly, as discussed in Section lTlI Bl However, the value 
of the degeneracy does not provide complete information 
about the density matrix of the source. In general, to 
characterize the source one needs to perform a complete 
state tomography. Without any assumptions about the 
source, this may become a practically impossible task. 



On the other hand, if the experimentator could make 
some reasonable assumptions about the source, based on 
his/her knowledge of the setup, then the state tomog- 
raphy might not be needed. For example, with an as- 
sumption about the structure of the density matrix, such 
as the one in Eq.®, one, in principle, can calculate the 
probabilities to detect vacuum, one-, and two-electron 
states (provided Nmax = 2) by fitting the plot of the to- 
tal degeneracy in Fig.(|4l). 

In certain situations, for instance, when the time res- 
olution of one's detector is not high enough, measuring 
the degeneracy as a function of time might not be practi- 
cal. Still, one would like to find probabilities of different 
components of the density matrix describing the source. 
Here is how it can be done with just one detector. 

First, the probability density to detect an electron, de- 
scribed by the state pi = | P^' ) ( 1 | in Eq. ([7]) , at a position 
X reads, 

p,{x,t) - {v^\¥{x,m{o\i>{x,t)\i^') 

where 4'(x,t) was previously defined in Eq.((22|) and, 

2 

as before, a(i) = m'^+^^AkH'^ ■ t = 

we immediately recover a familiar gaussian distribution 
Pi{x) = ■^^^ea:p(— 2AA:^x^) centered around .x = with 
the width l/2Afc. 

On the other hand, in the case of the two-electron state 
P2 in Eq. pH . the joint probability to detect one electron 
in the position xi a,t t = 0, and the other one in the posi- 
tion X2 a.i t = Q for symmetric/antisymmetric component 
is given by, 

Ps/As{xuX2) = |(2|'/^5|«'t(a;i)«'t(x2)|0)|2 

- .il±m?) " (l±-^^(^o-^o)), (32) 

where X = xi — X2- The total joint probability P(xi, X2) 
reads, 

1 3 

P{xi,X2) = -Ps{xi,X2) + -Pas{xi,X2). (33) 

The probability P{xi ,X2) ~ a two-detector quantity - can 
be used to construct a one detector probability P[x). 
For example, P{x) can be the probability to detect an 
electron(does not matter which) at a position x, pro- 
vided that the other electron is traced out. It can be 
obtained by integrating the probability P{xi,X2) over 
either one of the variables i.e. P{x) — dyP{x, y) — 

J^^dzP{z,x). Physically, P{x) can be measured with 
just one detector, by shutting it off once it has detected 
an electron(note that this is different from detecting both 
electrons at x; the probability for this event is given by 
P{x, x)). By integrating Eq. ([55)) over one of the variables 
we obtain, 

P(x) = ^Ps{x) + ^Pas{x), (34) 
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where the following quantities were introduced 

Ps/As{x) = (^|0)2) (^=^ ^^1'^) cosa;(fco - fco)). 

(35) 

Once again, Eq. represents the probability density to 
detect just one of the two electrons that are described by 
the state p2 while the destiny of the other electron is com- 
pletely ignored. If the overlap between single-electron 
states It/j) and \(f>) is non- negligible, then, by examining 
the Eg. (1351) . we conclude that P{x) will be reduced at the 
center position (x — 0) when compared to the probabil- 
ity Pi{x) in Ea ipri) . This is due to the interference terms 
containing — cosa;(fco — kg). The interference stems from 
the indistinguishability of electrons. When just one out 
of two electrons is detected, we cannot tell which electron 
was detected, and thus the interference occurs. 



P(x,y,t,t+T) 




y, (m) 

FIG. 5: Joint probability to detect one electron by a detector 
located at position {x, t) and another electron by a detector at 
(y, t + t) as a. function of the position of the second detector. 
The first detector is fixed at a; = 0.1m, t = Q, and t = 8.427 ■ 
10~®s. The parameters for the two-electron wave packet used 
in the simulations are ko^ — 1.024T0^^m~^, (fcoa, fco^) = for 
the first electron and k^^ — ka^ — 0.5- Afc, (fco^ , k^^ ) = for the 
second electron. The momentum spread Afc = 1.28 ■ 10*m~^ 
is the same for both electrons in all directions. 

The one-detector two-electron interference described 
here is a spatial version of the temporal behavior of 
the two-electron degeneracy component studied in Sec- 
tion HiTbI Both effects can be explained by the indis- 
tinguishability of individual electrons in a two-electron 
wave packet. These effects illustrate how two-electron 
phenomena are different from one-electron phenomena, 
and, thus, they can be naturally used to discriminate 
one- and two-electron contributions in the density ma- 
trix describing the source. 

Recall the density matrix ps describing the source and 
defined in Eq.Q. Imagine that from some experimental 



consideration we know that Nmax = 2, however, what is 
unknown are the probabilities of one- and two-electron 
outcomes pi,P2- Moreover, we only have one electron 
detector that we can move around such that a spatial 
distribution of electron detection events can be recorded. 
Of course, measuring the probability of having no elec- 
trons po is relatively easy since we just need to record 
no-click events and normalize their number to the total 
number of measurements. Thus, p2 can be expressed as 
a linear function of pi i.e. p2 = 1 — po — pi- Then, 
if we measure only one-electron events at different loca- 
tions and shut the detector off right away after it clicks, 
the resulting probability of one-electron events is given 
by piPi{x) + p2P{x), where Pi{x) and P{x) were de- 
fined in Eq. ([3T1) and Eq.lpl]) respectively. So, given the 
measured data, we can determine the value of pi and 
P2 by fitting the data with the one-parameter function 
piPi{x) + p2P{x), and, thus resolve the one- and two- 
electron components of the density matrix. 



B. Hanbury Brown and Twiss Effect for 
Two-electron States 

In many situations, one may be interested in coher- 
ence properties of the electron source (e.g. transverse 
(longitudinal) coherence, and (or) the effective size of the 
source) rather then in its density matrix. To determine 
those quantities, a two-detector measurement technique 
developed by Hanbury-Brown and Twiss (HBT) in stellar 
astronomy can be borrowed. 

In their groundbreaking work, HBT Q were interested 
in two-point intensity correlations of the light coming 
from a distant star. Using this data they were able to 
measure the angular diameter of the distant star. The 
joint two-detector probability P(xi, X2) given in Eq. ((33)) 
is the equivalent of the HBT measurement for pulsed 
electron sources at i = 0. Naturally, a pair of spa- 
tially separated electron detectors is required to measure 
this probability. Note that Eq. describes a particular 
measurement setup where one is interested in detecting 
two electrons in two different positions at the same in- 
stance in time(namely at t = 0). For electron sources 
whose electron statistics does not explicitly depend on 
time (e.g. continuous sources) Eg. (1551) also provide the 
probability of simultaneous detection of two electrons at 
any instance in time t > 0, provided that electron dis- 
persion effects are taken into account. For pulsed sources 
the expression in Eq. ([55)) has to be modified to account 
for time dependence in two-electron correlations at dif- 
ferent times (e.g. x,t and y^t + r). For instance, in 
Fig.®, the joint probability to detect one electron at 
{x — 0.1m, t = 8.427 • lO^^s) as a function of the posi- 
tion of the second detector at the same instance in time 
{t = 8.427 • 10~^s) is plotted. We observe that when the 
positions of both detectors coincide {x = y = 0.1m), the 
probability to find both electrons is reduced. This is a 
direct consequence of the PEP. The electrons that belong 
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to the antisymmetric two-electron state are not allowed 
to be in the same position. 



V. EFFECTS OF COULOMB REPULSION 

So far no interaction between electrons has been as- 
sumed. In reality, electrons in a two-electron state 
like, for example, the one described in Eq. p3|) are close 
enough spatially to experience strong Coulomb repulsion. 
But will this affect electron-electron interference and re- 
duce the degeneracy? Intuitively it is clear that the rela- 
tive momentum between two electrons will be increased 
by Coulomb repulsion, pushing them away from each 
other in momentum space. This, in turn, will result in a 
reduction of quantum interference. At first, one might, 
prematurely, conclude that all is lost and Coulomb re- 
pulsion would ruin quantum effects. However, if the rate 
of electron wave packet dispersion - a purely quantum 
feature - is larger than the rate of relative momentum 
change due to Coulomb repulsion then quantum effects 
will still dominate! 

Let us first estimate the effect of Coulomb interac- 
tion and then compare it to the rate of wave-packet dis- 
persion. For the former, consider two classical particles 
of equal charge and mass m,,. In the center-of-mass 
(COM) reference frame relevant system variables are the 
relative momentum p and the relative position of the par- 
ticles X. If the relative momentum of the particles at 
t = is po and the relative distance is xq then, by em- 
ploying the energy conservation argument, we arrive at 
the following relationship between x and p at t > 0, 



p 

m, 



(36) 



X me Xq 

which leads to the expression for p as a function of a;. 



P 



nieEo - 



mekeql 



(37) 



To proceed further, an estimate of the value of i?o is 
needed. We will assume that = h{ko^—k'Q^) = 0.5-ftAfc 
and we will use Afc = 1.28 • 10*m~^ - the same val- 
ues that we have used to describe two-electron states 
in the preceding section. The initial relative distance 
Xq can then be estimated from the distance between 
the maxima of the two-electron probability density as 
xo = -y^^ = 5.52 • lO^^m. On Fig.® the relative mo- 
mentum normalized to /i as a function of the relative 
distance starting from xq is plotted. As expected, the 
relative momentum increases with the relative distance 
until it reaches its maximum - a point where Coulomb 
repulsion is negligible. Using Eg. ([57)1 the terminal rel- 
ative momentum pt = y/m^EQ when a; — >■ cxd, and the 
ratio pt/po can both be determined. For our choice of 
parameters, this ratio is « 9. Next we calculate the 
relative distance Xd^t between the electrons when their 
COM reaches the detector placed O.lm away from the 
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FIG. 6: Relative momentum as a function of relative coor- 
dinate for two classical electrons (red solid line). The initial 
relative momentum is represented by the green solid line par- 
allel to the X axes. 



source by integrating Ea.([57|. After some algebra we 



find that Xdet — 5.77 • 10 



Note that the relative mo- 



mentum between electrons at Xdet is very close to its ter- 
minal value(see Fig.®). Finally, the effect of Coulomb 
repulsion for two electrons can be estimated by the ratio 
xa/xdet « 0.95 • 10-5. 

On the other hand, the rate of quantum dispersion of 
a two-electron wave packet in one dimension is given by 
the ratio Ax(0)/Ax(tdet) — \/ce{tdet), where a{t) was 
defined in Section IIII Al Calculating it for the same 
parameter values as in the Coulomb repulsion problem 
gives Ax{0) / Ax{tdet) ~ 3.1 • 10^^. Finally, by compar- 
ing magnitudes oixdet and Ax{tdet) we conclude that for 
our parameters Coulomb repulsion will, indeed, dominate 
quantum effects by roughly one order of magnitude. 

Fortunately, there is a way to enhance quantum effects 
while keeping Coulomb repulsion nearly constant. Note 
that the strength of Coulomb repulsion depends on the 
initial distance between electrons. The further away the 
electrons are, the smaller the effect of Coulomb interac- 
tion is. But simply putting the electrons further apart 
without increasing their longitudinal coherence Ax will 
again result in Coulomb domination. That means that to 
make quantum effects dominant it is necessary to increase 
electron coherence. Here is how that can be done in prac- 
tice. In reality, two-electron wave packets are three di- 
mensional. They are characterized by the longitudinal 
coherence Ax and transverse coherence Ay-Az. The lon- 
gitudinal coherence is related to coherence time At and 
is hard to control in an experiment. On the other hand, 
the transverse coherence can be easily improved by plac- 
ing a transversal filter in combination with a quadrupole 
lens. Increasing transverse coherence, while keeping lon- 
gitudinal coherence intact, will result in a wave packet 



11 



shaped like a pancake where longitudinahy the electrons 
are close, and trans versally far from each other. 

Consider for now a two-dimensional wave packet with 
the initial longitudinal coherence Aa; = 1/Afc, where 
Afc is the same as in the one-dimensional calculations. 
Assume that the initial transverse coherence Ay is ten 
times better than Ax, i.e. Ay = 10 • Ax. Like in 
the one-dimensional case, quantum effects can be esti- 
mated by the ratio Ax{0)Ay{0) / Ax{tdet)Ay{tdet) which 
is « 9.75 • 10~*. The effect of Coulomb repulsion can be 
also estimated by upgrading the one-dimensional model 
to two dimensions. The initial distance between the elec- 
trons is ro = 5.55-10~^m. The relative distance change at 
t — tdet due to Coulomb repulsion is r^/r^^^ « 9.35-10"^. 
And one immediately concludes that quantum effects are 
now of the same order as Coulomb repulsion. The quan- 
tum effects can further be made dominant by considering 
three-dimensional wave packets with improved transverse 
coherence. 



VI. SUMMARY AND CONCLUSIONS 

We developed a formalism that can be applied to the 
description of both pulsed and continuous nano tip elec- 
tron emission sources. To make quantitative predictions 
about particle statistics of a pulsed electron source, a 
model, describing the density matrix of the source, was 
introduced. Based on this model, the quantum degen- 
eracy of the source was calculated. The degeneracy at 



the source is predicted to be three orders of magnitude 
higher than the best previously reported value for contin- 
uous sources. Quantum effects due to electron-electron 
interference, that can be observed while measuring, for 
example, the degeneracy as a function of time were stud- 
ied. It is argued that these effects can readily be observed 
experimentally. The relationship to other electron inter- 
ference effects, such as the HBT effect, were discussed. 
The effect of Coulomb repulsion was estimated and it was 
shown that quantum effects can be made dominant over 
it. 

For the newest pulsed electron sources the expected 
vales of degeneracy are such that the Pauli Exclusion 
Principle should be considered to describe its properties 
correctly. Furthermore given the developments of tem- 
poral lenses to control and reverse electron dispersion, 
degeneracy is also expected to be essential at the image 
location of such temporal lenses. These developments 
appear to usher in a new field of electron quantum op- 
tics. Current efforts are underway in the laboratory [la] 
to observe such effects, and it is our hope that the above 
formalism provides some guidance to this and other sim- 
ilar future experiments. 
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